Computational Difficulty of Computing the Density of States 
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We study the computational difficulty of computing the ground state degeneracy and the density 
of states for local Hamiltonians. We show that the difficulty of both problems is exactly captured by 
a class which we call #BQP, which is the counting version of the quantum complexity class QMA. 
We show that #BQP is not harder than its classical counting counterpart #P, which in turn implies 
that computing the ground state degeneracy or the density of states for classical Hamiltonians is 
just as hard as it is for quantum Hamiltonians. 
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Understanding the physical properties of correlated 
quantum many-body systems is a problem of central im- 
portance in condensed matter physics. The density of 
states, defined as the number of energy eigenstates per 
energy interval, plays a particularly crucial role in this 
endeavor. It is a key ingredient when deriving many 
thermodynamic properties from microscopic models, in- 
cluding specific heat capacity, thermal conductivity, band 
structure, and (near the Fermi energy) most electronic 
properties of metals. Computing the density of states 
can be a daunting task however, as it in principle involves 
diagonalizing a Hamiltonian acting on an exponentially 
large space, though other more efficient approaches which 
might take advantage of the structure of a given problem 
are not a priori ruled out. 

In this Letter, we precisely quantify the difficulty of 
computing the density of states by using the powerful 
tools of quantum complexity theory. Quantum complex- 
ity aims at generalizing the well-established field of clas- 
sical complexity theory to assess the difficulty of tasks re- 
lated to quantum mechanical problems, concerning both 
the classical difficulty of simulating quantum systems as 
well as the fundamental limits to the power of quan- 
tum computers. In particular, quantum complexity the- 
ory has managed to explain the difficulty of computing 
ground state properties of quantum spin systems in vari- 
ous settings, such as two-dimensional (2D) lattices [l[ and 
even one-dimensional (ID) chains as well as fermionic 
systems Q. 

We will determine the computational difficulty of two 
problems: First, computing the density of states of a lo- 
cal Hamiltonian, and second, counting the ground state 
degeneracy of a local gapped Hamiltonian; in both cases, 
the result holds even if the Hamiltonian is restricted to 
act on a 2D lattice of qubits, or on a ID chain. To this 
end, we will introduce the quantum counting class ^BQP 
(sharp BQP), which constitutes the natural counting ver- 
sion of the class QMA (Quantum Merlin Arthur) which it- 
self captures the difficulty of computing the ground state 
energy of a local Hamiltonian [J H[ . Vaguely speaking, 
#BQP counts the number of possible "quantum solu- 
tions" to a quantum problem that can be verified using a 



quantum computer. We show that both problems, com- 
puting the density of states and counting the ground state 
degeneracy, are complete problems for the class #BQP, 
i.e., they are among the hardest problems in this class. 

Having quantified the difficulty of computing the den- 
sity of states and counting the number of ground states, 
we proceed to relate #BQP to known classical counting 
complexity classes, and show that #BQP equals #P (un- 
der weakly parsimonious reductions). Here, the complex- 
ity class #P counts the number of satisfying assignments 
to any efficiently computable boolean function. This can 
be a very hard problem which is believed to take exponen- 
tial time; in particular, it is at least as hard as deciding 
whether the function has at least one satisfying input, 
i.e., the complexity class NP. Examples for #P-complete 
problems (i.e., the hardest problems in that class) include 
counting the number of colorings of a graph, or com- 
puting the permanent of a matrix with binary entries. 
Phrased in terms of Hamiltonians, what we show is that 
computing the density of states and counting the ground 
state degeneracy of a classical spin system is just as hard 
as solving the same problem for a quantum Hamiltonian. 

Quantum complexity classes. — Let us start by intro- 
ducing the relevant complexity classes. The central role 
in the following is taken by the verifier V , which veri- 
fies "quantum solutions" (also called proofs) to a given 
problem. More formally, a verifier checking an n-qubit 
quantum proof (that is, a quantum state |^>)) consists of a 
T = poly(n) length quantum circuit U = Ut ■ ■ • U\ (with 
local gates Ut) acting on m = poly(n) qubits, which takes 
the n-qubit quantum state as an input, together 
with m — n initialized ancillas, \0)a = |0 • ■ ■ 0)a, applies 
U, and finally measures the first qubit in the {|0)i, |l)i} 
basis to return 1 ("proof accepted") or ("proof re- 
jected"). Then, the class QMA contains all problems of 
the form: "Decide whether there exists a such that 
Pucc{V(ip)) > a, or whether p acc (V(ip)) < b for all \yb), 
for some chosen a — b > l/poly(n), given that one is the 
case". Here, the acceptance probability of a state \tp) is 
P«cc(V(V0) := MJW. with 

n = (i/ <g> (ou)[/t(|i)(i|i ® i)u(i! ® \o) A ) , (i) 

which we illustrate in Fig. [TJ 
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FIG. 1: A QMA verifier consists of a sequence of T local uni- 
tary gates acting on the "quantum proof" \ip) and an ancillary 
register initialized to |0). The final measurement on the first 
qubit returns |1) or |0) to accept or reject the proof, respec- 
tively. Transition probabilities can be computed by doing a 
"path integral" over all intermediate configurations (ik)k- 



The idea behind this definition is that QMA quanti- 
fies the difficulty of computing the ground state energy 
Eq{H) of a local Hamiltonian H up to l/poly(n) accu- 
racy. Let the verifier be a circuit estimating (ip\H\ip); 
then a black box solving QMA problems can be used 
to compute Eq(H) up to l/poly(n) accuracy by binary 
search using a single QMA query. Note also that QMA is 
the quantum version of the class NP, where one is given 
an efficiently computable boolean function f(x) G {0, 1} 
and one must figure out if there is an x such that 

m = 1. 

The class NP has a natural counting version, known 
as ^P. Here, the task is to determine the number rather 
than the existence of satisfying inputs, i.e., to compute 
\{x : f(x) = 1}|. We will now analogously define #BQP, 
the counting version of QMA. Consider the verifying map 
51 of Eq. ([1]) for a QMA problem, with the additional 
promise that VL does not have eigenvalues between a and 
b, a — b > l/poly(n). Then the class #BQP consists of 
all problems of the form "compute the dimension of the 
space spanned by all eigenvectors with eigenvalues > a" . 

An equivalent definition for #BQP (cf. also [(| 0]) is 
the following: Consider a verifier D, with the additional 
promise that there exist subspaces A®1Z = C 2 such that 
(tp\ft\ip) > a for all \ip) e A, and (ip\n\i/>) < b for all \ip) € 
1Z, where again a — b > 1 /poly(n) - we can think of A and 
1Z as containing the good and bad witnesses, respectively. 
(Note that there will always be "mediocre" witnesses — 
the question is whether there exists a decomposition into 
a good and a bad witness space.) Then, #BQP consists of 
all problems of the form "compute dim A" . This number 
is well-defined, i.e., independent of the choice of A and 
1Z, and moreover, one can easily show that it is equivalent 
to the definition above, cf. the Supplementary Material. 

The gap promise we impose on the spectrum of fl is 
not present in the definition of QMA (though similarly 
restricted versions of QMA were defined in [y, 0]). Nev- 
ertheless, this promise emerges naturally when consider- 
ing the counting version: QMA captures the difficulty of 
determining the existence of an input state with accep- 



tance probability above a, up to a "grace interval" [b, a] 
in which mistakes are tolerated (i.e., if the largest eigen- 
value of f2 is in [6, a] , the oracle can return either out- 
come). Correspondingly, #BQP captures the difficulty 
of counting the number of eigenvalues above a, where 
eivenvalues in the grace interval [b, a] can be miscounted. 
The reason why we choose to define #BQP with a gap 
promise rather than with a grace interval is the same as 
for QMA, namely to have a unique outcome associated 
with any input. 

Similarly, the idea of the Hamiltonian formulation of 
the problem which we will discuss below is to ask for the 
number of eigenstates in a certain energy interval, where 
states which are in some small l/poly(n) neighborhood 
of this interval may be miscounted; again, for reasons of 
rigor we choose to consider only Hamiltonians with no 
eigenstates in that interval. It should be noted, however, 
that all of the equivalence proofs we give equally apply 
if we choose to allow for miscounting of states in those 
grace intervals instead of requiring them to be empty, as 
the proofs do not make use of the gap promise itself, but 
rather show that all states outside those grace intervals 
are mapped (and thus counted) correctly. Thus, while 
the actual number returned by the grace interval formu- 
lation of the counting problems might change under those 
mappings due to different treatment of states in the grace 
interval, it will still be in the correct range. 

The class #BQP inherits the important property from 
QMA of being stable under amplification, that is, the def- 
inition of #BQP is not sensitive to the choice of a and 
b. In particular, any a — b > l/poly(n) can be ampli- 
fied (by building a new poly-size fl' from Q) such that 
a' = 1 — exp(— poly(n)), b' = exp(— poly(n)), and keep- 
ing the eigenvalue gap between a' and by using a 
construction called strong amplification, cf. Rcf. Q; as 
shown there, strong amplification acts on all eigenval- 
ues independently and thus also applies to #BQP. The 
crucial point is that strong amplification works without 
changing the proof itself, compared to weak amplifica- 
tion which takes multiple copies of the proof as an input. 
While this is fine for QMA, it does change the dimension 
of the accepting subspace in an unpredictable way and is 
thus not an option for the amplification of #BQP. 

Complexity of computing the density of states. — Let 
us now show why the class #BQP is relevant for physical 
applications. In particular, we arc going to show that 
computing the density of states of a local n-spin Hamil- 
tonian H = ^jifj with few-body terms Hi, < 1, 
up to accuracy l/poly(n), is a problem which is complete 
for #BQP, i.e., it is as hard as any problem in #BQP can 
be. The same holds true for the (a priori weaker) prob- 
lem of counting the ground state degeneracy of a local 
Hamiltonian, given a l/poly(n) spectral gap above (note 
that Bravyi et al. suggested this as a definition for 
a quantum counting class). We can impose additional 
restrictions on the interaction structure of our Hamilto- 
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man, and as we will see, the hardness is preserved even 
for 2D lattices of qubits, or ID systems. 

The problem DOS (density of states) is defined as fol- 
lows: Given a local Hamiltonian H = compute 
the number of orthogonal cigenstates with eigenvalues in 
an interval [E\, E2] with E2 — E\ > l/poly(n), where we 
do not allow for eigenvalues within a small grace interval 
of width A = (E 2 — £'i)/poly(n) centerd around Ei and 
E 2 ; alternatively, we can allow for errorneous counts of 
eigenstates in that interval. Second, the problem #LH 
(sharp local Hamiltonian) corresponds to counting the 
number of ground states of a local Hamiltonian which 
has a spectral gap A = l/poly(n) above the ground state 
subspace, given we are told the ground state energy, and 
where we allow for a small splitting in the ground state 
energies; again, we can alternatively allow to miscount 
states in the grace interval. 

Clearly, #lh is a special instance of DOS, i.e., solving 
#lh can be reduced to solving DOS. In order to show 
that DOS is contained in #BQP, we can use a phase es- 
timation circuit (lpj to estimate the energy of any given 
input \ip) and only accept if its energy {ip\H\ip) is in the 
interval [E\, E2}) as the desired accuracy A = 1/poly (n), 
this can be done efficiently. A detailed proof (using a 
more elementary circuit) is given in the Supplementary 
Material. 

Let us now conversely show that #lh is a hard prob- 
lem for #BQP, that is, any problem in #BQP can be 
reduced to counting the ground states of some gapped 
local Hamiltonian [22J. As in turn #LH can be reduced 
to DOS, which is contained in #BQP, this proves that 
both #lh and DOS are complete problems for #BQP, 
i.e., they capture the full difficulty of this class. To this 
end, we need to start from an arbitrary verifier circuit 
U = Ut ■ ■ - U\ and construct a Hamiltonian which has as 
many ground states as the circuit has accepting inputs 
(corresponding to the outcome |l)i on the first qubit). 
Let A and 1Z be the eigenspaces of f2 [Eq. ([T])] with eigen- 
values > a = 1 - 2-P°W n ) and < b = 2-P ol y(™\ respec- 
tively, and define U[K] := {U\ip)i\0) A ■ € 71}. 

We will follow Kitaev's original construction for a 
Hamiltonian encoding a QMA verifier circuit [IHH, which 
for any proof \tp)i € A has the "proof history" |$) = 

2t=o ' ' ' Ui\ip}i\0}A\t)T as its ground state, where 
the third register is used as a clock. The Hamiltonian 
H = H init + J2j=i H cvoi{t) + fffinai has three types of 
terms: H init = 1 <g> (1 - |0)(0|a) <£> |0)(0|t makes sure the 
ancilla is initialized, H evo \(t) — — Ut <E) \t)(t — 1|t + h.c. 
ensures proper evolution from t — 1 to t, and -fffi, la i = 
^UVR] ® \T)(T\t gives an energy penalty to states |<fr) 
built from proofs \tp)i e 1Z. Note that our Hfi na i differs 
from the usual choice |0)(0|i ® 1 65 \T}(T\t and is in fact 
non-local; as we show in the Supplementary Material, 
this does not significantly change the relevant spectral 
properties (in particular, we keep the l/poly(n) gap, and 
the ground state subspace is split at most exponentially). 



With this choice of Hf Lna \ 1 H acts independently on the 
subspaccs spanned by {Ut • ■ ■ Ui\ip) i\x) A|i)r}t=o....,T for 
any \ip) € A or \ip) € 1Z, and \x)a the computational 
basis, and the restriction of H to any of these subspaces 
describes a random walk which is characterized by the 
choice of \tp)i and the number of l's in |aj)>i. These cases 
can be analyzed independently (see Supplementary Ma- 
terial) , and it follows that H has a dim .4-fold degenerate 
ground state space with a 1/poly (rt) gap above, proving 
#BQP-hardness of #lh. 

This shows that #lh is #BQP-hard for a Hamiltonian 
which is a sum of log T- local terms (i.e., each term acts 
on log T sites), as the clock register is of size log T. In 
order to obtain a fc-body Hamiltonian, Kitaev suggested 
to use a unary encoding of the clock (i.e., \t)r is encoded 
as 1 1 - - - 10 - ■ - 0) , with t l's), so that each Hamiltonian 
term only acts on three qubits of the clock. However, 
this makes the Hilbert space of the clock too big, and 
terms need to be added to the Hamiltonian to penal- 
ize illegal clock configurations. These terms divide the 
Hilbert space into two parts, "Hi eg ai©'Hiiicgai- Here, W\ cg &\ 
contains only legal clock states, whereas HiiiRgaiContains 
only configurations with illegal clock states [4|, |5[ . Since 
no Hamiltonian term couples these two subspaces, the 
Hamiltonian can be analyzed independently on the two 
subspaces. It turns out that its restriction to Hiiicgai has 
an at least l/poly(n) higher energy, while on 7Li em \, the 
Hamiltonian is exactly the same as before. Thus, one 
finds that the new Hamiltonian still has the right num- 
ber of ground states, and a l/poly(n) spectral gap. The 
very same argument applies in the case of ID Hamiltoni- 
ans, using the QMA-construction of Ref. Again, the 
Hamiltonian acts independently on a "legal" and an "ille- 
gal" subspace, where the latter has a polynomially larger 
energy, and the former reproduces the (low-energy) spec- 
trum of the original Hamiltonian (llj . 

An alternative way to prove QMA- hardness on re- 
stricted interaction graphs is to use so-called perturbation 
gadgets, which yield the Hamiltonian of the Kitaev con- 
struction above from a pcrturbative expansion; in partic- 
ular, this way one can show QM A-hardness of Pauli-type 
nearest-neighbor Hamiltonians on a 2D square lattice of 
qubits [H- As shown in Ref. 12 1, such gadgets do in fact 
approximately preserve the whole low-energy part of the 
spectrum, and thus, our #BQP-hardness proof for #LH 
still applies to these classes of Hamiltonians. It should 
be noted, however, that since the eigenvalues are only 
preserved up to an error l/poly(n), the splitting of the 
ground state space will now be of order 1/poly (n); how- 
ever, it can still be chosen to be polynomially smaller 
than the spectral gap. 

Quantum vs. classical counting complexity. — As we 
have seen, the quantum counting class #BQP exactly 
captures the difficulty of counting the degeneracy of 
ground states and computing the density of states of local 
quantum Hamiltonians. In the following, we will relate 
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#BQP to classical counting classes and prove that #BQP 
is equal to #P, counting; the number of satisfying inputs 



to a boolean function [23l | . In physical terms, this shows 



that counting the number of ground states or determining 
the density of states for a quantum Hamiltonian is not 
harder than either problem is for a classical Hamiltonian. 

Clearly, #P is contained in #BQP, as the latter in- 
cludes classical verifiers. It remains to be shown that 
any #BQP problem can be solved by computing a #P 
function. We start from a verifier operator fl, Eq. ([T]), 
and wish to show that the dimension of its accepting 
subspace, i.e., the subspace A with eigenvalues > a, can 
be computed by counting satisfying inputs to some ef- 
ficiently computable boolean function. Using amplifica- 
tion, we can ensure that |dim.A — tr fi| < |, i.e., we need 
to compute tr O to accuracy j . This can be done using 
a "path integral" method, which has been used previ- 
ously to show containments of quantum classes in the 
classical classes PP and #P (see e.g. [HI). We rewrite 
tr O = /(C) as a sum over products of transition prob- 
abilities along a path C = (io, . . . , ijv, ji, • • • , Jjv), where 

/(C) =<io|/<OU^|ji)0'il^ ■ • • 4lJT>x (2) 
(i T \ [\0)(0\i ® 1] |«tX*tI^t ' ■ • Ui\i )l\0) A 

(cf . Fig. [IJ . Since any such sum over an efficiently com- 
putable /(C) can be determined by counting the satis- 
fying inputs to some boolean formula (see the Supple- 
mentary Material for details), it follows that SI can be 
computed using a single query to a black box solving #P 
problems. 

Summary and discussion. — In this work, we consid- 
ered two problems: Computing the density of states and 
computing the ground state degeneracy of a local Hamil- 
tonian of a spin system. In order to capture the com- 
putational difficulty of these problems we introduced the 
quantum complexity class #BQP, the counting version 
of the class QMA. We proved that this complexity class 
exactly captures the difficulty of our two problems, even 
when restricting to local Hamiltonians on 2D lattices of 
qubits or to ID chains, since all these problems are com- 
plete problems for the class #BQP 24 1. 

We have further shown that #BQP is no harder than 
its classical counterpart #P. In particular this implies 
that computing the density of states is no harder for 
quantum Hamiltonians than it is for classical ones. While 
this quantum-classical equivalence might seem surpris- 
ing at the Hamiltonian level, it should be noted that 
the classes #P and PP quite often form natural "upper 
bounds" for many quantum and classical problems. 

What about the problem of computing the density of 
states for fcrmionic systems, such as many-electron sys- 
tems? On the one hand, this problem will be still in 
#BQP and thus #P, since any local fermionic Hamilto- 
nian can be mapped via the Jordan- Wigner transform to 
a (non-local) Hamiltonian on a spin system, whose energy 



can still be estimated efficiently by a quantum circuit [3] ■ 
On the other hand, hardness of the problem for #BQP 
can be shown e.g. by using the #BQP-hardness of #lh, 
and encoding each spin using one fcrmion in two modes, 
similar to 14|. Thus, computing the density of states for 
fermionic systems is a #BQP-complete problem as well. 
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Note added. — After completion of this work, we 
learned that Shi and Zhang [15] have independently de- 
fined #BQP and shown its relation to #P using the same 
technique. 



[1 
[2 

[3 

[j: 

[5 
[6 

[7; 

[8 
[9 
[10 

[11 

[12 

[13 
[14 
[15 



[16 
[17 

[18 
[19 



R. Oliveira and B. M. Terhal, Quant. Inf. Comput. 8, 



900 (2009), quant-ph/0504050 



D. Aharonov, D. Gottesman, S. Irani, and J. Kempe, 
Commun. Math. Phys. 287, 41 (2009). larXiv:0705.4077l 
N. Schuch and F. Verstraete, Nature Physics 5, 732 
(20091. larXiv:0712.0483l 

A. Y. Kitaev, A. H. Shen, and M. N. Vyalyi, Classical and 
quantum computation (American Mathematical Society, 
Providence, Rhode Island, 2002). 

D. Aharonov and T. Naveh, (2002), [quant-ph/02100"77l 
D. Aharonov, M. Ben-Or, F. G. Brandao, and O. Sattath, 
(2008) , larXiv:0810.4840l 

R. Jain, I. Keren idis, M. Santha, and S. Zhang, (2009), 
larXiv:0906.4425l 

C. Marriott and J. Watrous, Comput. Complex. 14, 122 
(2005), [cs70506068| 

S. Bravyi, C. Moore, and A. Russell, (2009), 
larXiv:0907.T297l 

R. Cleve, A. Ekert, C. Macchiavello, and M. Mosca, Proc. 



R. Soc. Lond. 454, 339 (1998), [^a^ph/9708016 



N. Schuch, I. Cirac, and F. Verstraete, Phys. Rev. Lett. 
100, 250501 (20081. [arXiv:0802.335ll 
J. Kempe, A. Kitaev, and O. Regev, SIAM Journal of 
Computing 35, 1070 (2006), |quant-ph/0406180l 
L. Adleman, J. DeMarrais, and M.-D. Huang, SIAM J. 
Comput. 26, 1524 (1997). 

Y.-K. Liu, M. Christandl, and F. Verstraete, Phys. Rev. 
Lett. 98, 110503 (2007), [o^iant-ph/0609125| 
Y. Shi and S. Zhang, private communication; 
notes available at http://www.cse.cuhk.edu.hk/ 
~syzhang/papers/SharpBQP .pdf . 
R. Bhatia, Matrix Analysis (Springer, 1997). 
C. Jordan, Bulletin de la Societe Mathematique de 
France 3, 103 (1875). 

P. Halmos, Trans. Amer. Math. Soc. 144, 381 (1969). 
E. Bernstein and U. Vazirani, Siam J. of Comp. 26, 1411 

(1997), [o3Iant-ph/970100l| 

[20] F.G.S.L. Brandao, PhD thesis (20081. larXiv:0810.0026l 



5 



[21] T. J. Osborne, (2006), |«md-mat/0605194"| 
[22] Note that the connection between QMA with a unique 
"good witness", such as in our #BQP definition, and 
local Hamiltonians with unique ground state and a 
l/poly(n) gap has been shown in [y]. 
[23] Formally speaking, the reduction from #BQP to #P is 
weakly parsimonious, i.e., for any function / 6 #BQP 
there exist polynomial-time computable functions a and 



P, and a function g £ #P, such that / = a o g o j3. This 
differs from Karp reductions where no postprocessing is 
allowed, a = Id, but it still only requires a single call to 
a #P oracle. 

[24] Note that computing the density of states to multiplica- 
tive accuracy is less difficult, see Refs. [2p|. |2lj. 



SUPPLEMENTARY MATERIAL 



Quantum Complexity Classes 

We start with the relevant definitions. Let x be a bi- 
nary string. Then, wc denote by the verifier V = V x a 
quantum circuit of length T = poly(|x|), U = Ut • • -U\ 
(with local gates Ut) acting on m = poly(|a;|) qubits, 
which is generated uniformly from x. The verifier takes 
an n = poly(|a;|) qubit quantum state \4>}i as an in- 
put (we will express everything in terms of n instead 
of | a; | in the following), together with m — n initialized 
ancillas, \0)a = |0---0)a, applies U, and finally mea- 
sures the first qubit in the {|0)i,|l)i} basis to return 
1 ("proof accepted") or ("proof rejected"). The ac- 
ceptance probability for a proof \ip) is then given by 
PaccTO)) := (-0|O|V>, with 

n = (lj ® <0U)C/+(|TXl|i ® 1)17(1/ ® |0) A ) • (3) 

Definition 1 (QMA). Let < a,b < 1 s.th. a - b > ^ 
for some polynomial p(n) . A language L is in QMA(a, b) 
if there exists a verifier s. th. 

xeL => A 

max 

x £ L X max (Cl) < b . 

Definition 2 (#BQP). Let < a, b < 1 s.th. a -b > 

l/poly(n), and let fl be a verifier map with no eigenvalues 
between a and b. Then, the class ^BQP(a, b) consists of 
all problems of the form "compute the dimension of the 
space spanned by all eigenvectors of H with eigenvalues 
> a" . 

An alternative definition for #BQP is the following: 

Definition 3 (#BQP, alternate definition). Consider a 
verifier il with the property that there exist subspaces A® 
TZ = C N (N = 2 n ) such that (V>|ft|^) > a for all \if>) G A, 
and (ip\n\ip) < b for all |^) € TZ, where again a — b > 
l/poly(n). Then #BQP(a, b) consists of all problems of 
the form "compute dirndl". 

Note that dim A is well-defined: Consider two decom- 
positions C N = A © TZ and C N = A' ® TZ' . Without loss 



of generality, if we assume dim .4 > dim A', it follows 
dim„4 + dim 7?.' > N, and thus there exists a non-trivial 
G AnTZ', which contradicts the definition. 

Theorem 4. Definition [H and Definition [3] are equiva- 
lent. 

Proof. To show that Definition[2]implies Definition[3l let 

A be spanned by the eigenvectors with eigenvalues > a. 
To show the converse, we use the minimax principle for 
eigenvalues [l6j], which states that the /cth largest eigen- 
value Afe of a Hermitian operator f2 in an iV-dimcnsional 
Hilbcrt space can be obtained from either of the equiva- 
lent optimizations 

Afc(Q) = max min (arlfila;} (4) 

Mk \x)£M k 

= min max (a;|f2|x) , (5) 

M N -k + i \x)£M N -k + i 

where Mk is a subspace of dimension k, and \x) is a unit 
vector. Now notice that Def. [3] implies that 

min (x|f2|x) > a and max (x|f2|x) < b. (6) 

\x)eA \x)en 

Next, consider the minimax theorem for k = dirndl. 
From Eq. Q we have 

Adim^t = max min (x|f2|x) 

> min (x|n|x) 
\x)eA 

> a. 

Now consider the case that k = dim„4+ 1. From Eq. (|5|), 
using the fact that ./V — (dim .4 + 1) + 1 = dim TZ, we have 

Adim.A+i = m i n max (x|f2|a;) 

< max(:r|f2|x) 

\x)EK 

< b. 

Thus we have 

XdimA > a > b> Adim.A+1 , (J) 
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since a — b > l/poly(n). It follows that Adim^t is the 
smallest eigenvalue of fi which is still larger than a, and 
therefore the span of the first dim A eigenvectors of fl 
is equal to the span of all eigenvectors with eigenvalue 
> a. The equivalence follows. □ 

The class #BQP(a, b) inherits the useful property of 
strong error reduction from QMA: the thresholds (a, 6) 
can be amplified to (1 — 2~ r ,2~ r ), r = poly(n) while 
keeping the size of the proof: 

Theorem 5. Let #BQP„(a,&) denote #BQP with an 
n qubit witness. Then ^BQP rl (a, b) C #BQP n (l — 
2~ r ,2~ r ) for every r € poly(n). 

This follows directly from the strong amplification pro- 
cedure presented in [8(, which describes a procedure to 
amplify any verifier map such that any eigenvalue 
above a (below b) is shifted above 1 — 2~ r (below 2~ r ) at 
an overhead polynomial in r. 

Using Thm. [SJ we will write #BQP instead of 
#BQP(a, b) from now on, where a — b > poly(n), and 
a, b can be exponentially close to 1 and 0, respectively. 



E 2 . Then, the problem ff=LH = ff=LH(Ei — Eq) (sharp 
local Hamiltonian) is to compute the dimension of the 
eigenspace with eigenvalues < E\. 

Note that #lh depends on the "energy splitting" 
E\ — Eq of the low-energy subspace. In particular, for 
E\ — Eq = 0, ^lh(O) corresponds to computing the de- 
generacy of the ground state subspace. As we will sec 
in what follows, the class #lh(ct) is the same for any 
splitting cxp(— poly(n)) < a < l/poly(n). 

We now show that #lh and DOS are both #BQP- 
complete. We do so by giving reductions from 
#LH(l/poly(n)) to DOS, from DOS to #BQP, and from 
#BQP to #Lti(cxp(— poly(n))) ; this will at the same 
time prove the claimed independence of #lh(ct) of the 
splitting a. 

Theorem 9. #LH(l/poly(n)) reduces to DOS. 

Proof. If we denote the parameters of the #lh problem 
by Eq , Ei , E2 , then we can simply relate them to the 
parameters E\, E 2l A of a DOS problem by A = E 2 — E \ , 
Ei = Eq — iA and E 2 = E\ + \ A, and the result follows 
directly. □ 



The Complexity of the Density of States 



Theorem 10. DOS is contained in #BQP. 



We now use the class #BQP to characterize the com- 
plexity of the density of states problem and the problem 
of counting the number of ground states of a local Hamil- 
tonian. We start by defining these problems, as well as 
the notion of local Hamiltonian, and then show that both 
problems are #BQP-complete. 

Definition 6 (k- local Hamiltonian). Given a set of 
poly(n) quantum spins each with dimension bounded by 
a constant, a Hamiltonian H for the system is said to 
be k-local if H = £\ Hi is a sum of at most poly(n) 
Hermitian operators Hi, \\Hi\\ < 1, each of which acts 
nontrivially on at most k spins. 

Note that fc-local does not imply any geometric local- 
ity, only that each spin interacts with at most k — 1 other 
spins for any given interaction term. However, we re- 
strict ourselves to k = 0(log(n)) so that each Hi can be 
specified by an efficient classical description. 

Definition 7 (dos). Let E 2 - £f > l/poly(n), A = 
(E%— £7i)/poly(n), and let H = J^. Hi be a k-local Hamil- 
tonian such that H has no eigenvalues in the intervals 
[Ei — A , Ei + A] an d [E 2 -§,E 2 + %]. Then, the prob- 
lem DOS (density of states) is to compute the number 
of orthogonal eigenstates with eigenvalues in the interval 
[Ei,E 2 ]. 

Definition 8 (#lh). Let Eq < E x < E 2 , withE 2 -E x > 
l/poly(n), and let H = Hi be a k-local Hamiltonian 
s.th. H > Eq, and H has no eigenvalues between E\ and 



Proof. We start with a fc-local Hamiltonian H as in 
Dcf. [7J Now define a new Hamiltonian 

H' := v{H 2 - (Ei + E 2 )H + E X E 2 ) . 

W is a 2fc-local Hamiltonian; here, v = l/poly(n) is 
chosen such that each term in H' is subnormalizcd. Any 
eigenvalue of H in the interval [Ei + A ; E 2 — A] translates 
into an eigenvalue of H' which is below 

A := -v±(E 2 - Ei - A) < -l/poly(n) , 



A . 



whereas any eigenvalue outside \E\ 
lates into an eigenvalues of H' above 



Eo 



trans- 



B :=u^{E 2 - Ei 



A ■ 
2 



> l/poly(n) . 



The original DOS problem now translates into counting 
the number of eigenstates of H' with negative energy, 
given a spectral gap in a I/poly (n) sized interval [A; B] 
around zero. We now use the circuit which was used in [5( 
to prove that log-local Hamiltonian is in QMA; it accepts 
any input state l^) with probability 

i (m'w 

pirn) = o ' 

where m = poly(?i) is the number of terms in H' . [The 
idea is to randomly pick one term H[ in the Hamiltonian 
and toss a coin with probability (1 — (ip\Hl\tp))/2.] 
Computing the answer to the original DOS problem 
then translates to counting the number of states with 
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acceptance probability > a = 5+2^, given that there 
are no eigenstates between a and b = \ + and since 
a — b = (A — B)/2m > l/poly(n), this shows that this 
number can be computed in ^BQP. □ 

Theorem 11. #LH(exp(— poly(n))) is #BQP-/iard 

Proof. To show #BQP-hardncss of #lh, we need 
to start with an arbitrary QMA verifier circuit U = 
Ut ■ ■ ■ U\ and construct a Hamiltonian with as many 
ground states as the circuit has accepting inputs. By 
amplification, we can assume that the acceptance and re- 
jection thresholds for the verifier are a = 1 — e and b = e, 
where we can choose e = 0(exp(— p(n))) for any polyno- 
mial p(n) . As before, let A and 72 be the eigenspaces of 
57 with eigenvalues > a and < 6, respectively. Define 

U[R] := {U\i/>)i\0)a : H')i € 72} (8) 

and denote the projector onto this space by Ilm^y . No- 
tice that for any state |x) G U[7l), due to our rejection 
threshold b — e, we have 

( X |(|l)(l|i®l)|x)<e. (9) 

We now follow Kitaev's original construction to encode 
a QMA verifier circuit into a Hamiltonian which has a 
"proof history" as its ground state for any proof \(f>)i G 
A [1, That is, the ground states of the Hamiltonian 
are given by 

T 

\9) = ^2u t ...Ui\4>)i\o) A \t) T (10) 

t=o 

where the third register is used as a "clock" . The Hamil- 
tonian has the form 

T 

H = H init + 2^ Hevol(t) + -fffinal (H) 
t=l 

where 

• Hinit = 1/ ® (1 - |0)(0|a) ® |0)(0| T checks that 
the ancilla is property initialized, penalizing states 
without properly initialized ancillas; 

• flevol(t) = -\U t <Z>\t){t-l\T-\Ul®\t-l)(t\ T 

+5I ® \t)(t\ T + 5I (8 \t- l)(t - 1| T 

checks that the propagation from time t — 1 to t is 
correct, penalizing states with erroneous propaga- 
tion; 

• -Hfinai = Hu[r.] ® \T}(T\t causes each state \4>) built 
from an input 1 G 72. (but which nonetheless has 
a correctly initialized ancilla) to receive an energy 
penalty. 



As we show in Lemma [12l the total Hamiltonian H has 
a spectral gap l/poly(n) above the dim ^4-dimensional 
ground state subspace. However, -ff nna i is not a local 
Hamiltonian, but as we argue in the following, it can be 
replaced by the usual term = \0)(0\i <8> 1 ® \T)(T\ T 

while keeping the ground space dimension (up to small 
splitting in energies) and the l/poly(n) spectral gap. As 
we prove in Lemma [T3l 

Hlnil > tffinal " ■ (12) 

Thus, replacing 7Jfi na i by Hf^ al will decrease the energy 
of any excited state by at most y/e = 0(exp(—p(n)/2)). 
(The energy of the ground states is already minimal and 
cannot decrease.) On the other hand, the energy of any 
proper proof history for H cannot increase by more than 

{x\HgUx) = (Xl|0)(0|i ® l|x) < e = 0(exp(-p(n)) , 

due to our choice of acceptance threshold, i.e., the low 
energy subspace has dimension dim A We thus obtain 
a Hamiltonian with a dim^ dimensional ground state 
subspace with energy splitting < e = exp(— p(n)) and a 
l/poly(n) spectral gap above. □ 

The following two lemmas are used in the preceding 
proof of Theorem QT] 

Lemma 12. H has a spectral gap of size l/poly(n). 

Proof. Our proof follows closely the discussion in 
Ref. [3], cf. also [H|]. We can block diagonalize H by 
the (conjugate) action of the following unitary operator, 

T 

W = 52Uj---Ui®\j){3\ T , (13) 
o=o 

which maps H — > W = W^HW. As this has no effect 
on the spectrum, we can work with the simpler H' from 
now on. 

Let us explicitly write the effect of conjugation by W 
on the terms of H. The first term is unaffected, H[ nit = 
TJinit. The final term becomes -ff^ nal = H^ <g> |0)(0|yi ® 
\T)(T\t, where H^. is simply the projector onto the space 
72, and the ancillas are in the correct initial state. 

We can now conjugate each of the terms in H evo i(t) 
separately. For example, the first term gives 

W\U t g> \t)(t - 1\)W = 1 ® !*><* — 1| . (14) 

The other terms are exactly analogous, and we find that 

H' eyol {t) = l®\[\t-l){t-l\ + \t){t\ 

— !*><* — II — I* — 1><*|] - 

The total evolution Hamiltonian is then block-diagonal 
a matrix which looks like a hopping Hamiltonian in the 
clock register, 

5^1^(0 = (is) 

t 
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where the (T + l)-by-(T + 1) -dimensional tri-diagonal 
matrix E is given by 



Proof. We begin by__ recalling the result due to Jor- 
dan [I?} (see Ref. 



18j for a more modern treatment) 



E 



\ 



1 - 



V 



(16) 



for the simultaneous canonical form of two projectors. 
In the subspace where P and Q commute, both opera- 
tors are diagonal in a common basis and the spectrum 
is either (0, 0), (0, 1), (1, 0), or (1,1), and direct sums of 
those terms. In the subspace where they don't commute, 
the problem decomposes into a direct sum of two-by-two 
blocks given by 



To discuss the spectrum of H', let 

51 =A®\0)(0\ A ®C T+1 , 

5 2 =11®\0)(0\ A <Z>C T+1 , 

and l>3 the orthogonal complement of Si © S2', i-e., Si 
corresponds to evolutions starting from good proofs, 1S2 
to those starting from wrong proofs, and 1S3 to evolutions 
with wrongly initialized ancillas. Note that H 1 = H\s 1 © 
H\s 2 (B H\s 3 , thus, we can analyze the spectrum for H\s p 
separately. Note that the restriction to S p does not affect 
H' evol (t). Since we expect the ground state subspace to 
occur on <Si, we need to compute ground state energy 
and gap of H'\ Sl , as well as lower bound the ground 
state energies of H'\g 2 and H' s . 

First, if/ nit |s = -HfinJs = 0, i.e., the spectrum of H'\ s 
equals the spectrum of E, which can be straightforwardly 
determined to be 1 — co; 



Pi 



1 





Qj 



c c 
cs s 



(18) 



c 2 + s 2 = 1, and the 



where s = sin(0j) for some angle 
subscript j just labels a generic block. 

In fact, this two-by-two block form is completely gen- 
eral if we allow embedding our projectors into a larger 
space while preserving their rank. The rank-preserving 
condition guarantees that our bound is unchanged, since 
we are only appending blocks of zeros, and so we will 
consider this two-by-two form without loss of generality 

The constraint that \\Q(t — P)Q\\ 00 < e implies con- 
straints on the values that sm(9j) can take. In particular, 
we can directly compute this operator norm in each block 
separately, and we find that for all j 



\Q 3 (1-P 3 )Q 3 \ 



\0 3 ) < 



cigenfunctions (cos i< 



with u n — mr/ (T + 1), and We can also directly compute in each block that 



); the ground state 



degeneracy is indeed dim A as desired. 

On the other hand, -fffi nal |s 2 = 1, and H[ nit \ Sa > 1, i.e., 
in both cases the ground state energy of H'\g p is lower 
bounded by the ground state energy of 



P j -Q j 



l-d- 1- 

—cs — s 



(19) 



(20) 



E' 



' 2 



\ 



V 



which has eigenvalues 1 — cosi9„, with i?„ = (n+i)7r/(T+ 
|), and eigenfunctions (cos \"& n , cos |t?„, . . . ). 

It follows that H' (and thus H) has a ground 
state energy of 1 — coswo = 0, and a gap 



1 

□ 



cos 



2T+3 



0(1/T 2 ) 



(9(l/poly(n)) above. 



It remains to prove Eq. (|12[) . which follows from the 



following Lemma by choosing P = |0) (0|i®l, Q = ILjpj] , 
and using Eq. (j9|). 

Lemma 13. Let P and Q be projectors such that ||Q(1 — 
P)Q\\co < e- Then 



P-Q> -y/el. 



(17) 



The spectrum of this operator is easily computed to 
be ±|sin(0j-)|. Thus, the least eigenvalue of P — Q is 
bounded from below by —y/e. □ 



Quantum vs. Classical Counting Complexity 

We finally want to relate #BQP to the classical count- 
ing class #P. It is clear that any #P problem can be 
solved in #BQP by using a classical circuit. We will 
now show that conversely, #BQP can be reduced to #P 
under weakly parsimonious reductions: That is, for any 
function / e ^BQP there exist polynomial-time com- 
putable functions a and fi, and a function g G #P, such 
that / = a o g o p. This differs from Karp reductions 
where no postprocessing is allowed, a = Id, but still only 
requires a single call to a #P oracle, in contrast to Turing 
reductions. 

Theorem 14. There exists a weakly parsimonious re- 
duction from #BQP to #P. 

Proof. We start from a verifier operator f2 for a #BQP 
problem. First, we use strong error reduction to let a = 
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1 - 2-("- 2 ) and b = 2- {n+2 \ It follows that 



Idim^-trfil < 2"2-(" +2) 



(21) 



and thus we need to compute tr f2 to accuracy 1/4. This 
can be done using the "path integral" method previously 
used to show containments of quantum classes inside PP 
and #P [13|| . We rewrite tr Q = 5If /(C)> where the sum 
is over products of transition probabilities along a path, 
which we label 



C = {io, ■ ■ ■ ,in, ji, ■ ■ ■ Jn) , 



(22) 



so that 



/(C) =(zo|/(OUC/ 1 t | J1 )(j 1 |C/ 1 t • • • 4|i T )x 

(i T | [10X01, <g> 1] \i T )(i T \U T ■ ■ ■ U^iohflU 



(23) 



(cf. Fig. 1 in the main manuscript for an illustration). 

Since any quantum circuit can be recast in terms of real 
gates at the cost of doubling the number of qubits 
we can simplify the proof by assuming /(C) to be real. 
To achieve the desired accuracy it is sufficient to approx- 
imate / up to |C| + 2 digits, where |C| = poly(n) is the 
number of bits in C- Now define 



3 (C) :=round[2^+ 2 (/(C) + l)] 



(24) 



and note that g(() is a positive and integer- valued func- 
tion satisfying 



[2- |c| - 2 £s(c)-i]-E/(o 



< 



(25) 



Finally, by defining a boolean indicator function, 



MOO 



i ifo<c<. 9 (C) 

otherwise 



we can write g(C) = X^>o MOO; an d thus, 

£s(C) = £M00 ■ 

c u 



This shows that tr Q can be approximated to accuracy 
j, and thus dim A can be determined by counting the 
number of satisfying assignments of a single boolean 
function /i(00 that can be efficiently constructed from 
f2, i.e., by a single query to a black box solving #P 
problems, together with the efficient postprocessing 
described by Eg. (f25|) . □ 



